Robust d x 2_ y 2 pairing symmetry in high-temperature superconductors 
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Although initially quite controversial, it has been widely accepted that the Cooper pairs in op- 
timally doped cuprate superconductors have predominantly d x i_ y i wavefunction symmetry. The 
controversy has now shifted to whether the high-T c pairing symmetry changes away from opti- 
mal doping. Here we present phase-sensitive tricrystal experiments on three cuprate systems: 
Yo.rCao.sBasCuaOr-s (Ca-doped Y-123), La 2 -*Sr x Cu0 4 (La-214) and Bi 2 Sr 2 CaCu208+5 (Bi- 
2212), with doping levels covering the underdoped, optimal and overdoped regions. Our work 
implies that time-reversal invariant, predominantly d x 2_ y 2 pairing symmetry is robust over a large 
variation in doping, and underscores the important role of on-site Coulomb repulsion in the making 
of high-temperature superconductivity. 



In the cuprate perovskites, doping with electrons 
or holes converts antiferromagnetic insulators into 
high-temperature superconductors with a normal-to- 
superconducting transition temperature T c that depends 
sensitively on the amount of doping Q. The doping de- 
pendence of T c can be described for many hole-doped 
cuprate systems by the empirical formula 0: 



T c (p) = T c , max [l-82.6(p-p c )% 



(1) 



where p is the number of holes per Cu02 layer and p c « 
0.16. Optimal doping represents a watershed in many su- 
perconducting and normal state properties of the cuprate 
superconductors 1]. Nanometer-scale charge inhomo- 
geneity, pseudogap, and other anomalous normal- state 
properties are observed mostly in the underdoped re- 
gion, while a Fermi-liquid description appears valid in 
the overdoped region jll. Moreover, the Hall number Q, 
superfluid pair density |4j and many other properties also 
exhibit remarkable dependencies on doping. In the lan- 
guage of quantum criticality |5|, there are a number of 
competing states near the quantum critical point (~ p c ). 
Variation in doping may induce a symmetry breaking in 
favor of, e.g. a spin density or charge density wave phase, 
or to a superconducting phase with different pairing sym- 
metry such as the time reversal symmetry broken pair 
states d x 2_ y 2 + is or d x 2_ y 2 + id xy |6|. 

Although d x 2_ y 2 pairing symmetry is well estab- 
lished for several optimally doped cuprate supercon- 
ductors |3, there are a number of indirect symme- 
try studies suggesting a doping-induced change in pair- 
ing symmetry in some cuprates. For example, tun- 
nelling spectroscopy suggests a significant gapped (s- 
wave) component in the pairing wavefunction in over- 
doped Yi_ ;E Ca :c Ba2Cu307_5 8], and a change in sym- 
metry from d x 2_ y 2 to d x 2_ y 2 + id xy or d x 2_ y 2 + is in 
overdoped YBa2Cus07 One penetration depth mea- 
surement in electron-doped Pr2-^Ce iC Cu04 indicated a 



d to s transition in the optimal to overdoped range 
jl0l | . although a second study indicated d pairing at 
all doping levels in this system ^]|. In addition, the 
low-temperature limit of the thermal conductivity ex- 
pected for a d-wave superconductor Il2l was observed 
in La2-^Sr ;E Cu04 for all doping levels 

Phase sensitive tests of pairing symmetry using the 
tricrystal geometry have been described elsewhere Q. 
Briefly, thin films of cuprate superconductors are epitax- 
ially grown on a SrTiOs substrate composed of 3 grains. 
The substrate and subsequent cuprate thin film geometry 
are chosen (Fig. 1(a)) such that, for a d x 2_ y 2 supercon- 
ductor, in the absence of supercurrents there are an odd 
number of sign changes in the component of the pairing 
wavefunction normal to the grain boundaries upon cir- 
cling the tricrystal point. The sign changes at the grain 
boundaries cost Josephson coupling energy. This energy 
is reduced by the generation of circulating supercurrents, 
resulting in a Josephson vortex with a vortex number 
= 1/2. For real superconducting order parameters and 
the conventional I c =Iisin(0) Josephson current-phase re- 
lationship the tricrystal vortex should have exactly half 
of the conventional flux quantum </>o=h/2e=2.07x 10 -15 
Wb of magnetic flux. In our geometry (Fig. 1(a)), a 
spontaneous =1/2 vortex will occur at the tricrystal 
point only if the d x 2_ y 2 component is at least as large as 
a possible s component. Further, a fixed imaginary com- 
ponent to the order parameter would make the absolute 
value of the flux in the tricrystal point vortex different 
upon inversion by roughly the ratio of the imaginary com- 
ponent to the total pairing amplitude. The =±1/2 
vortices at the tricrystal point could have the same abso- 
lute values of flux, if there were domains with imaginary 
order parameter components, but such domains would 
have to alternate signs on a length scale smaller than 
the experimental spatial resolution. In our experiments, 
the flux at the tricrystal point is imaged with a scanning 



Superconducting Quantum Interference Device (SQUID) 
microscope (SSM). In our SSM, the sample is scanned 
using a mechanical lever mechanism relative to a well 
shielded pickup loop integrated into a SQUID sensor |l4j . 

Fig. l(c,d,e) shows SSM images of =±1/2 vortices 
at the tricrystal point of a 130nm thick Ca-doped Y-123 
film epitaxially grown by pulsed laser deposition on a 
SrTiOs tricrystal in the geometry of Fig. la, chosen to 
show the half-flux quantum effect for a superconductor 
with d x 2_ y 2 pairing symmetry. These films have suffi- 
ciently high grain boundary supercurrent densities [0| 
that Josephson vortices in the grain boundaries and at 
the tricrystal point are resolution limited with a 4 mi- 
cron diameter pickup loop. The =+1/2 state can be 
inverted to its degenerate time-reversed -1/2 state (Fig. 
1(c)) by applying currents of a few mA with the appro- 
priate polarity through the field coil of the SQUID sus- 
ceptometer [lq . 

Integration of the total flux at the tricrystal point (Fig. 
1(f)) shows that the = ±1/2 states have the same to- 
tal flux and field distribution. The time reversal symme- 
try broken pair states such as d x 2_ y 2 J r id xy and d x 2_ y 2 J r is 
are thus ruled out. 

In contrast to Ca-doped Y-123, La-214 films grown on 
SrTiOs tricrystals have relatively low T c 's and very small 
grain boundary critical current densities, presumably due 
to the lattice mismatch between the deposited film and 
the substrate 17]: The fields of the half-fluxon at the 
tricrystal point for the T C =28.5K sample shown in Fig. 
2 spread out along the grain boundaries over tens of mi- 
crons. Nevertheless, it is clear that there are spontaneous 
currents at the tricrystal point, and that the resulting 
vortex retains time- reversal symmetry (Fig. 2(b)). For 
detailed modelling, the flux/unit length in the ith branch 
of the vortex can be written as (18) 

dn 2tt Aji 1 + a?e" 2r */ A i< 1 ) 

where n is the absolute value of the distance along 
the ith grain boundary from the tricrystal point, J, 
is the Josephson penetration depth of the ith grain 
boundary, ^/(Aj^l + af)) is the same for each grain 
boundary, and the total flux at the tricrystal point 
is (j) = Xli(^o/27r)4 tan _1 (ai). The flux through the 
SQUID pickup loop is calculated by propagating the 
fields at the surface to the height of the pickup loop using 
Fourier transform techniques, and then integrating over 
the known pickup loop geometry. Results of fits of these 
calculations to the data (Fig. 2(c)) show that to within 
our experimental errors, the vortex at the tricrystal point 
has a total flux of 4>q/2. 

Tricrystal experiments were also performed on a num- 
ber of Bi-2212 superconductors with various doping lev- 
els achieved by controlling the oxygen content during film 
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FIG. 1: (a) Tricrystal geometry. The polar plots represent 
the pairing wavefunctions, with the red lobe phases shifted 
by 7r relative to the blue. (b) Schematic of the pickup 
loop area of the SQUID susceptometer sensor used. (c,d,e) 
are SQUID microscope images of the tricrystal point of a 
Yo.7Cao.3Ba2Cu307-5 film on a SrTi03 tricrystal with the 
geometry of (a). All SSM images in this paper were taken 
at 4.2K. The dashed lines indicate grain boundaries, (c) il- 
lustrates the inversion of a =-1/2 vortex at the tricrystal 
point (top image), (c, middle image) a 5mA pulse of current 
is passed through the susceptometer field coil to invert a - 
1/2 vortex (c, top image) to a ±1/2 vortex, creating also a 
= -1 Josephson vortex in the horizontal grain boundary, 
(c, bottom image) the -1 Josephson vortex is dragged from 
the tricrystal point by moving the sensor parallel to the grain 
boundary while applying a current of 4mA. The color scales 
span <j) s — 0.66 0o of flux through the SQUID pickup loop in 
(c) and 0.38 0o in (d) and (e). (f) Shows the integration of 
the total flux (in units of <fio) of the = ±1/2 state (d, red 
dots), the -1/2 state (e, blue dots), and a nearby = +1 
integer vortex (green dots) over a circular area Ai n t centered 
at the tricrystal point. The blue line in (f) is the =-1/2 
data multiplied by -1, demonstrating time reversal invariance. 
The red line in (f) is the = ±1/2 data multiplied by 2. 




(b) 




(c) ToO/Im 



1 

• Data 


1 


1 


_ — Model 
- B 




B' - 




•v \ 

7 \. 




A 

*^ sss!SfSS ^ 1 


1 


A' 



-100 -50 50 100 

Position (/im) 

Aj,AA' = 36/LUT1 A J,BB' = 30 /^ m 



FIG. 2: SQUID microscope images, taken with an 8/im 
square pickup loop, of the tricrystal point region for an under- 
doped (T C =28.5K, thickness=310nm) La-214 film epitaxially 
grown by laser ablation, (a) Difference image (0 s (±2OnT)- 
s (-2OnT))/2 (to cancel out stray fields) for images taken with 
the sample cooled in ±20nT fields, resulting in N0 = ±1/2 
Josephson vortices at the tricrystal point. The sum image 
(0 s (+2OnT)+0 s (-2OnT))/2 (b) shows that the flux from the 
= ±1/2 vortices cancel out within a few percent, demon- 
strating time-reversal symmetry, (c) Cross-sections through 
the data of (a) along the lines indicated by the arrows, and 
fits using the Josephson penetration depths Aj along the two 
grain boundaries as fitting parameters, assuming the vortex 
at the tricrystal point has 0o/2 of magnetic flux. The best fit 
value, allowing the total flux at the tricrystal point to vary, 
was 0=O.585±O.l fa. 



deposition, or thermal annealing of a given tricrystal Bi- 
2212 film (Fig. 3). Modelling like that described 
above [HI showed that in all cases the amount of mag- 
netic flux at the tricrystal point was 0o/2 to within exper- 
imental error, and that the absolute value of the total flux 
at the tricrystal point remained the same within a few 
percent when the vortex was inverted, either by applying 
a local field at low temperature, or by cooling again in 
a slightly different field. The crystal structure of all the 
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FIG. 3: SQUID microscope images for tricrystal samples with 
optimally, underdoped, and overdoped Bi2Sr2CaCu20s+5 
300nm thick films epitaxially grown by RF sputtering on a 
SrTiOs substrate with the geometry of Fig. 1(a). The images 
were taken with a 4/zm octagonal pickup loop (optimal) and a 
7.5/xm square pickup loop (overdoped and underdoped). Full 
scale variation in <j) s was O.40o, O.180o, and 0.21 0o respec- 
tively. The lines are cross-sections of the image data through 
the tricrystal point along the directions indicated. In each 
case there is a = +1/2 vortex at the tricrystal point. 

Bi-2212 samples studied here is tetragonal equivalent Q. 
A d+s mixed pair state in Bi-2212 is therefore symmetry 
forbidden. 

Fig. 4 displays the doping range covered in this study. 
Numerous studies indicate that the crystal structure in 
multi-grain cuprates is preserved except in a narrow re- 
gion lnm in width along the grain boundaries |2£|. It 
is possible that the dopant concentration at the grain 
boundary is different from the bulk. However, our 
temperature dependent scanning susceptibility measure- 
ments on Yo.7Cao.3Ba2Cu307_ ( 5 bicrystals indicate that 
the supercurrent across the grain boundary has a T c 
within 5K of that in the bulk, and it is unlikely that 
the dopant profile near the grain boundary remains con- 
stant while the grains themselves undergo a full range 
of variation. Therefore, our results show that in several 
cuprate systems, over a wide range of doping, pairing 
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YBa2Cii307_5, suggesting that the d £L 2_ y 2 nodal quasi- 
particles survive to very low doping |3(l |. 

The present work, coupled with previous work estab- 
lishing (i-wave symmetry for a number of optimally hole- 
and electron-doped cuprates Q , calls for a universal ori- 
gin of d x 2_ y 2 pairing symmetry in all the cuprate super- 
conductors studied so far. We suggest that all evidence 
points to a strong influence of strong on-site Coulomb 
repulsion, a characteristic common to all cuprates, and 
which is also responsible for the doping-induced Mott 
insulator-metal transition observed in all cuprate per- 
ovskites. 
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FIG. 4: Plots of T c /T c?max vs. p, the calculated dop- 
ing holes per CuC>2 layer, using Eq. 1, with p c = 0.16, 
and T c ,max = 90K, 90K, 38K and 82K for optimally doped 
YBa2Cu307-<5 (yellow dot), Yo.7Cao.3Ba2Cu307-<5 (green 
dot), La2-ccSr x Cu04 (red dots) and Bi2Sr2CaCu20s+(5 (blue 
dots) respectively. 

with d x 2_ y 2 symmetry is robust, and that any fixed imag- 
inary component of the order parameter in pair states 
such as d x 2_ y 2 + id xy and d x 2_ y 2 + is must be small. 
Note that d x 2_ y 2 pairing symmetry in the Bi-2212 sys- 
tem persists to p ^0.07 (Fig. 4), very close to the onset of 
superconductivity at p=0.05. At such low doping other 
orders such as anti-ferromagnetism, charge density wave 
and spin density wave phases compete vigorously with 
d-wave superconductivity 0, E| • 

As Anderson |2l|, |22| first suggested, the physics 
of both the normal and superconducting states of the 
cuprates is dictated by strong Coulomb interactions 
in the Cu02 planes, although the exact form of the 
ground-state wavefunction is still a matter of debate |2^ . 
Many theoretical studies including strong correlation ef- 
fects favour d x 2_ y 2 over d xy and extended s- wave (s*) 
states H |2g, with pairing in the d x 2_ y 2 channel en- 
hanced by on-Cu-site repulsion and suppressed by inter- 
site Coulomb interactions 26]. The stability of the pure 
d x 2_ y 2 pair state is further enhanced by the presence of 
the van Hove singularity or flat band around (Q-ZLh i^fi) 
in the 2D band structure of the Cu02 planes @,|24[. The 
origin of such Fermi surface pinning can be attributed to 
the effect of on-site repulsion |27j |. 

The robust nature of d x 2_ y 2 pairing over a wide doping 
range (0 < p < 0.35) has been demonstrated by several 
numerical studies based on Hubbard models [23, |221, con ~ 
sistent with our present work. Furthermore, our observa- 
tion of d x 2_ y 2 pairing in the low doping regime p w 0.07 
is supported by a recent measurement of c-axis penetra- 
tion depth as a function of temperature and doping in 
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